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We theoretically calculate the finite wave vector plasmon dispersion in a low density 2D electron layer 
taking into account finite temperature, finite layer width, and local field corrections. We compare our 
theoretical results with recent Raman scattering spectroscopic experimental 2D plasmon dispersion 
data in GaAs quantum wells at very low carrier densities (rs > 10) and large wave vectors (g > 
kp)- We find good agreement with the experimental data, providing an explanation for why the 
experimentally measured dispersion seems to obey the simple classical long wavelength 2D plasmon 
dispersion formula. We also provide a critical discussion on the observable manifestations of the 
quantum-classical and the Wigner crystal - electron liquid crossover behavior in the 2D plasmon 
properties as a function of electron density and temperature in GaAs quantum well systems. 
PACS Number : 73.20.Mf, 73.21.-b, 71.45. Gm 



A plasmon (or, a plasma excitation) is a fundamental 
elementary excitation of an electron system (or, for that 
matter, any charged particle system). It is the collec- 
tive (normal) mode of charge density oscillation in the 
free carrier system, which is present both in classical and 
quantum plasmas. Studying the collective plasmon ex- 
citation in the electron gas has been among the very 
first theoretical quantum mechanical many body prob- 
lems studied -Ln solid state physics dating back to the 
early 1950sErcl. A great deal of theoretical and exper- 
imental work has been carried out during the last fifty 
years on the issue of observable many body effects in 
the plasmon dispersjon relation, originally in the con- 
text of bulk metalscl ajid more recently for two dimen- 
sional (2D) plasmonsd occurring in 2D electron gases 
confined in semiconductor inversion layers, heterojunc- 
tions, quantum wells, and superlattices. The specific 
theoretical issue is that the plasmon frequency (wp) is 
exactly known only at long wavelengths {q — ^ 0) where 
the /-sum rule fixes the plasma frequency to be neces- 
sarily that given by classical electrodynamics, namely, 
LOp — (27rne^/m)^/^Y/g in 2D and ujp = (47rne^/m)^/^ 
in 3D where n is the 2D (3D) electron density and 
is interpreted as / n where k is the background 

lattice dielectric constant. At finite g, away from the 
long wavelength limit, there are several corrections to 
the plasmon dispersion ujp{q) arising from nonlocal (finite 
wave vector) response, finite temperature thermal cor- 
rections, many-body effects, local field corrections, and 
other mechanisms relevant to the specific electron system 
being studied (e.g., finite width of the 2D layer, band 
structure effects, interface effects, mode coupling effects 
to phonon and other possible elementary excitations, im- 
purity scattering or disorder effects, etc.). Significant 
deviations from tbsplong wavelength classicai, result are 
predicted to occurErt2l in 2D (as well as in 3Dcl — in this 
paper we restrict our considerations, motivated by recent 
experimentalj, entirely to 2D plasmons) plasmon disper- 
sion as the electron density is lowered, and the plasmon 



wave vector ceases to be in the long wavelength limit 
{q <C kp, where kp = {Ann/gy^^, with g as the degen- 
eracy factor of the 2D electron system including spin, is 
the 2D Fermi wave vector) and starts becoming compa- 
rable {q « kp) to the Fermi wave vector. Using the usual 
dimensionless electron gas (quantum) coupling strength 
parameter (= I/qb^/t^, where as {= fi^/rne^) is the 
effective Bohr radius) , where rs is the ratio of the nonin- 
teracting kinetic energy to the average potential energy 
at T = 0, one expects large finite wave vector correc- 
tions to the classical long wavelength plasmon disper- 
sion relation ujci{q) — (^o^/q, with = 2TTne'^/m, for 
^ 1 where the electron gas is strongly interacting. In 
fact, there are theoretical predictions, based^on model 
many-body plasmon dispersion calculationstj, that the 
2D plasmon frequency may renormalize to zero at finite 
wave vectors (and rg 3> 1) due to interaction-induced 
softening of the plasmon dispersion. In general, theories 
including higher order many-body and other effects pre- 
dict a lowering of the plasma frequency below the classi- 
cal plasmon dispersion curve for rg > 1 and finite values 
of q/kp. 

It therefore comes as a surprise that recent inelastin 
light (Raman) scattering based direct measurementstiJ 
of 2D plasmon dispersion in high quality low density 
(down to electron densities as low as n « 5 x lO^cm"^ 
corresponding to rg parameter as high as 25) electron 
systems in GaAs quantum wells find remarkably good 
agreement between the measured 2D plasmon disper- 
sion and the simple classical formula {ujp « lUci) upto 
wave vectors as Ipge as q ^ 2kp. Earlier experimen- 
tal measurementstS of 2D plasmon dispersion in semi- 
conductor structures were typically restricted to higher 
carrier densities (rs < 1) and lower wave vectors {q ^ 
0.1 — 0.2kp) where the classical q^^^ plasmon disper- 
sion relation is valid by virtue of the applicability of the 
long wavelength /-sum rule. A primary motivation for 
the recent loai, electron density 2D plasmon dispersion 
measurementtj, which we are trying to theoretically un- 
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derstand in this paper, has been to explore high-wave 
vector (q > kp) and low density (rg 3> 1) dispersion cor- 
rections to the classical q^/"^ 2D_plasmon dispersion for- 
mula. The unexpected findingtll that the experimental 
2D plasmon dispersion follows quantitatively the classi- 
cal q^/"^ formula upto large values of g/fci? (~ 2) even for 
very strongly interacting (r^ ~ 25) 2D electron systems is 
a significant puzzle, particulaelv^a view of the extensive 
existing theoretical literatureOil^l on finite wave vector 
many-body and nonlocal plasmon dispersion corrections 
showing very large deviations from the classical plasma 
frequency. 

By carrying out a realistic random phase approxima- 
tion (RPA) calculation of the finite temperature, finite 
wave vector 2D plasmon dispersion in the actual GaAs 
quantum wells used in ref. we provide a partial res- 
olution of the puzzle posed by the data presented in ref. 
[ri| . In particular, we show that RPA provides an excel- 
lent quantitative description of the data — the nonlocal 
finite q corrections introduced by RPA (which are sub- 
stantial for q/kp ~ 2) are not explicitly manifest in the 
experimentally measured plasmon dispersiontll because 
of a fortuitous cancellation among a number of indepen- 
dent contributions to the plasmon energy, mostly the al- 
most exact accidental numerical calculation between the 
higher-order (i.e., g'^/^ and higher) nonlocal finite wave 
vector correction (tending to increase the plasma fre- 
quency above the classical dispersion curve and becom- 
ing quantitatively prominent for q/kp > 1) and the finite 
layer width correction arising from the quasi-2D nature of 
the electron layer (which by itself lowers the plasma fre- 
quency below the classical dispersion curve because the 
usual q^^ Coulomb interaction in a strictly 2D system is 
softened in the quasi-2D system, and this softening in- 
duced lowering of the plasma frequency is particularly 
quantitatively significant in low density systems, where 
the layer width in the transverse direction is large due to 
the weakening of the Hartree selfconsistent potential con- 
tribution to the confinement). We also find, as explained 
below, that a slight 10%) adjustment in the quoted 
sample densities in ref. |ll|, where density has not been 
independently measured but inferred from the measured 
plasmon dispersion, significantly improves the agreement 
between our theoretical results and experimental data. 

The plasmon dispersion can be calculated by finding 
poles of the density-density correlatipn-iunction. Within 
RPA (or its simple generalizationsQclll3 including local 
field corrections arising from correlation effects) plasmon 
modes at finite wave vectors and finite temperatures are 
given by the zeros of the dielectric function, e(<7,a;;T) = 
1 — f ((7)no(<7, w; T), where v{q) is the Coulomb interac- 
tion modified by both the quasi-2D form factor due to 
the finite width of the 2D quantum well and correlation- 
induced local field effects, and llo{q,u!,T) is the nonin- 
teracting 2D finite temperature irreducible polarizability. 
The modification of the Coulomb interaction due to cor- 
relation induced local field correction is modeled by a 
static correlation factor G{q), which we calculate within 



the Hubbard approximationpOM. In the presence of local 
field corrections the Coulomb interaction is modified in 
the following manner: v{q) ^'('?)[1 ~ G{q)]. Note that 
local field corrections, < G{q) < 1, tend to soften the 
Coulomb interaction because exchange-correlation effects 
tend to keep the electrons away from each other reducing 
the effective Coulomb interaction. We calculate numeri- 
cally the plasmon dispersion by solving e(g, lo;T) — to 
obtain u)p{q] n, T) by incorporating thermal, finite thick- 
ness, and local field correlation effects. (All other effects, 
e.g. phonons, give negligible corrections to the plasmon 
dispersion in a low density GaAs quantum well electron 
system.) For each density the chemical potential of the 
2D system has to be calculated self-consistently at tem- 
perature T since at the low densities of interest to us, 
where Ep kgT, the chemical potential is very differ- 
ent from the Fermi energy. In our calculation we incor- 
porate finite thickness corrections by using a quasi-2D 
form factor appropriate for an infinite square well poten- 
tial with a width d. At the low electron densities used 
in ref. |ll| the infinite square well confinement model is a 
very good approximation since the Hartree selfconsisijcpt. 
effects are weak. We include the local-field correctionQ't2l 
using the Hubbard approximation (HA), where the 2D 
HA correlation factor is given by 

where ko = kp at T = 0. We generalize the T = HaS-EhI 
to finite temperatures by simply redefining ko to be the fi- 
nite temperature analogy of the Fermi wave vector using 

the formula: ko{T) = kp{T/Tp)ln{l + e^/'^^^), where 
^ = ^(T) = A:BTln[e("/'=^^)(27./gm) _ ^ ^j^^ gj^j^^, 

perature 2D chemical potential obtained from the total 
number of particles n — J dED{E)f{E), using the 2D 
electron density of states, D{E), and the Fermi distribu- 
tion function, f{E). The use of this particular (i.e., HA) 
local field correction is an uncontrolled approximation of 
our theory since there is no generally accepted theoreti- 
cal framework to incorporate correlation effects (beyond 
RPA) in the plasmon dispersion calculation. Our finite 
temperature generation of the HA is, however, quite rea- 
sonablft-since it correctly interpolates between the T = 
2D HAQ with ko{T — 0) = kp and the high temperature 
result, ko(T/Tp — > oo) = 0, of vanishing local field cor- 
rection. We have checked that the use of other staticJpcal 
field corrections (e.g, the finite-T STLS correctionst5) |dp| 
not significantly alter our conclusion. It is well-knowru'Ll 
that the HA for the local field correction, while being 
static and theoretically uncontrolled, has the great ad- 
vantage of being simple and reasonably quantitatively 
accurate compared with other local field correction ap- 
proximations. 

In Fig. 1 we show our calculated plasmon dispersion 
along with the available experimental data from ref. |l^. 
In Figs. 1(a) and 1(b) respectively plasmon dispersion 
curves for rg — 5.2 {n = 1.2 x 10^"cm^^) and = 17.2 
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(n = 1.1 X lO^CTO"^) are shown at T = 200mK. We 
use in our calculations system parameters corresponding 
to GaAs-Ala;Gai_2.As quantum wells with a well width 
d = 330A as appropriate for ref. |l|. Thin solid lines rep- 
resent the classical long wavelength plasmon dispersion 
(lUci = ujQy^). Dotted lines show the dispersion includ- 
ing finite wave vector non-local (i.e., higher order in q) 
effects calculated within RPA, which leads to an increase 
in plasma frequency compared with ujd ■ The leading or- 
der RPA dispersion correction is given by ujp{q)/ujci — 
1 + (3/4)(g/gTF), where qtf (= gme^/h^) is the 2D 
screening wave vector. Dashed lines are the dispersions 
including corrections by both non-local and finite thick- 
ness effects. Finite well width reduces the plasma fre- 
quency by softening the 2D Coulomb interaction, and its 
long wavelength dispersion in an infinite square well po- 
tential model can be calculated by straightforward alge- 
bra to be given by ujp{q)/ujci = 1 — 0.207{d/aB){q/qTF)- 
Noting that g = 2 (spin degeneracy) in GaAs quantum 
wells, i.e. qxpOLB = 2, we find that the two dispersion cor- 
rections (non-local and finite width) approximately can- 
cel each other when iq/Aqxp = 0-207dq/aBqTF- Thus 
we find that non-local effects and finite thickness cor- 
rections cancel each other almost exactly in a quan- 
tum well with well thickness d « 3. 605. The GaAs 
(ttB ~ 90 A) sample used in ref. |Tl| has a quantum well 
width d = 330A, which gives rise to an almost exact 
cancellation between non-local effects and finite thick- 
ness corrections. As shown in Fig. 1 RPA results with 
realistic finite thickness effects agree with the classical 
formula upto a large wave vector [q ^ 2kp) and provide 
reasonable quantitative agreement with experimental 2D 
plasmon dispersion data. Dot-dashed lines in Fig. 1 are 
the plasmon dispersion curves including local field cor- 
rections, non-local dispersion, and finite thickness effects 
at r = 0. The local field corrections tend to reduce 
the plasma frequency since it softens the Coulomb in- 
teraction by keeping the electrons effectively away from 
each other, and in the long wavelength limit the 2D plas- 
mon dispersion within the HA local field corrections is 
given at T = by LUp{q)/uJci = 1 - (rs/2V2)(g/(?TF)- 
As the density decreases (or rg increases) the reduction 
of the plasma frequency due to local field corrections is 
stronger, and this effect dominates in large wave vectors. 
Thick solid lines show the full plasmon dispersion includ- 
ing finite temperature effects at T = 2QQmK including all 
the corrections described above. (Long dashed line repre- 
sents the plasmon dispersion including all corrections at 
T = 200mK with zero temperature HA, i.e., /cq = ^f-) 
Even within RPA finite temperature effect by itself in- 
creases the plasma frequency. The thermal enhancement 
is negligible (exponentially small for T/Tp ^ 1) at low 
T in high density GaAs 2D electron systems since the 
Fermi temperature in these systems is much greater than 
the experimental temperature {T/Tp <C 1). However, a 
very low density electron system (e.g., electrons on the 
surface of liquid Helium, where the Fermi temperature is 
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FIG. 1. Calculated plasmon dispersions with available 
experimental data [11] for (a) n = 1.2 x 10^°cm~^ [r^ = 5.3) 
and (b) n = 1.1 x 10** cm"^ (r^ = 17.2). Thin solid line is the 
classical local plasmon dispersion and thick solid line indicates 
the dispersions with all corrections (as explained in the text) 
at a finite temperature (T = 2QQmK). Other lines in the 
figure are described in the text. The arrows indicate the Fermi 
wave vector and Fermi energy. The dots are experimental 
data points from ref. [11]. 

smaller than or comparable with experimental temper- 
ature) is essentially a non-degenerate classical plasma, 
where the long wavelength plasmon dispersion is given 
by ujp{q)/ujci = 1 + {3/2\/2){T/Tp){q/rs). In a quantum 
plasma (T <C Tp) on the other hand the leading order 
thermal correction to the T = plasmon dispersion is 
exponentially weak in temperature. We thus see that for 
T /Tp not too small (which is the experimental situation 
in ref. |ll|, where T/Tp m 0.4 at the lowest density) the 
thermal correction to the plasma frequency tends to can- 
cel the local field correction similar to the cancellation 
between the nonlocal finite q corrections and quasi-2D 
finite width corrections. 

In Fig. 1(a) the plasmon dispersion for rg = 5.3 
(n — 1.2 X 10^°cm~^) is shown. In this relatively high 
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density system the experimental data lie in the long wave- 
length {q < kp = 2.7 X lO^CTO"^) and the low temper- 
ature (T ^ Tp = 5.0K) limit. Note that the enhance- 
ment of plasma frequency by non-local effects is almost 
cancelled by finite thickness effects. The local field cor- 
rections and finite temperature effects are not important 
in this sample since these effects are quantitatively signif- 
icant only at low n or large rg. In Fig. 1(b) we show the 
plasmon dispersion for — 17.2 (n = 1.1 x lO^cm"^). 
Experiment datallJ for this sample lie in the effective large 
wave vector {q > kp ~ 0.83 x lO^cm^^) regime. Even in 
this low density sample (and at large effective wave vec- 
tors) the experimental data can be seen to agree very well 
with the classical plasma dispersion and the finite thick- 
ness RPA calculations. Since the non-local effects are 
almost canceled by the finite width correction, we spec- 
ulate that the reduction of the plasma frequency due to 
local field corrections is perhaps almost exactly canceled 
by the thermal enhancement of the plasma frequency. 
Our calculated results (Fig. 1(b)), however, show the lo- 
cal field corrections to be too large in this low density 
sample to cancel out with thermal effects. It is certainly 
possible, perhaps even likely, that our calculated HA local 
field corrections are strong overestimations of the actual 
finite temperature [T/Tp ~ 0.4 in the experiment) local 
field corrections, and in reality local field corrections, be- 
ing much smaller that our HA results, do in fact cancel 
out with the finite temperature enhancement. This re- 
mains an important open question for future theoretical 
work. One would need a quantitatively more accurate 
and theoretically well-controlled finite temperature the- 
ory for local field corrections for this purpose. At present 
no such local field theory exists in the theoretical litera- 
ture. 

Since the local field corrections (even within HA) are 
regarded as improvement to RPA in calculating many 
physical properties of low density systems (strongly cor- 
related systems) the large discrepancy between the plas- 
mon mode dispersion (in Fig. 1(b)) including local field 
effects and the experimental data in this low density elec- 
tron system is unexpected. It is troublesome to uncriti- 
cally accept that correlation effects are negligible based 
on the mere speculation that finite temperature local field 
corrections may be small (in other words, much smaller 
than what we calculate within the HA). We now suggest 
another possibility. In ref. |ll] the electron densities are 
estimated using the classical local plasmon dispersion for- 
mula based on the fact that the ^ dispersion seems to 
apply very well to the experimental plasma dispersion. 
This could, however, be problematic, and may in fact 
lead to an underestimation of n. We now compare the 
zero temperature classical local plasmon formula used to 
estimate the electron densities given in ref. |ll| (rt = 1.1 
and 3.7 x lO^cm"^) with the corresponding finite tem- 
perature plasmon dispersion including all the corrections 
described above, using somewhat higher densities in the 
compete dispersion calculation. In Fig. 2 we show the 
plasmon dispersion including all corrections for densities 
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FIG. 2. Comparison of finite temperature (T = 2QQmK) 
plasmon dispersions (solid lines) including all corrections 
for the density n = 1.37,4.1 x lO" cm~^ with zero 
temperature classical local dispersions (dashed lines) for 
n = 1.1,3.7 x 10^ cm~^. Experimental data are taken from 
ref. [11]. Inset shows the plasma frequency as a function of 
density parameter at g = A;f. Dot dashed line is the clas- 
sical plasmon mode, dashed lines are RPA modes (thick line 
for T — 300m_ft' and thin line for T = OA'), and solid lines 
are HA modes (thick line for T = 300mK and thin line for 
T = OK). 

n = 1.37 (lower solid line) and 4.1 x lO^cm"^ (upper 
solid line) . Experimental data points are taken from ref. 
pi] . Dashed lines in Fig. 2 are the classical T — local 
plasmon dispersion for the density n = 1.1 (lower line) 
and 3.7 x lO^cm^^ (upper line). Our calculated full plas- 
mon dispersion in Fig. 2 for somewhat higher densities 
including all corrections also agrees very well with ex- 
periment just as the classical T = formula apparently 
does for the lower densities proposed in the experiment. 
Thus, whether the local field corrections are large or not 
is not obviously clear until one can measure the exper- 
imental electron density using an independent method. 
This is, however, very difficult to do in the interesting 
regime of very low electron densities where the usual Hall 
density measurement fails. At this stage all we can say 
is that the full plasmon dispersion (including local field 
corrections and the other effects) could explain the exper- 
imental data provided one uses somewhat larger electron 
densities than the experimental carrier density estimates 
(based on a comparison of the data with the classical 
T — plasma dispersion formula) in ref. |ll|. Inset in Fig. 
2 shows the quantum-classical crossover of the plasma 
frequency at wave vector q = kp and T = 300mK as a 
function of the density parameter rg . In the high density 
limit the enhancement of the plasma frequency by non- 
local effects dominates all other corrections, but in the 
low density limit the local field corrections give rise to a 
large decrease of the plasma frequency. Thermal correc- 
tion of the plasmon mode also increases as the density 
decreases. 
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FIG. 3. Comparison of plasmon dispersions calculated 
within HA (solid lines, the thick and thin Une being for 
T = 200mK and T = OK, respectively) with the phonon 
modes of a hexagonal Wigner lattice (dashed lines) for 
Ts = 25. The dotted line is the classical plasmon dispersion. 

In Fig. 3 we compare the plasmon dispersion calcu- 
lated within the quantum HA with the "phonon" modes 
of the electron solid (hexagonal Wigner crystal (WC)o) 
within the harmonic approximation upto the zone bound- 
ary of the lattice for — 25 (corresponding to an elec- 
tron density, n ~ 5.4 x lO^cm^^, for the 2D GaAs sys- 
tem). Inset shows the hexagonal Wigner lattice in the 
momentum space, lul (i^t) indicates the longitudinal 
(transverse) phonon mode of the 2D Wigner lattice. The 
longitudinal WC phonon mode (which corresponds to the 
plasmon in the electron liquid system) has much lower 
frequency than the quantum plasmon mode at = 25, 
which is consistent with the fact that the transition to 
a WC phase is expected to occur at substantially lower 
electron densities (around rg = 37 or below). In the long 
wave length limit the WC phonon mode agrees with plas- 
mon mode within HA. But in the high wave vector region 
q> kp {kp = 5.7 x lO^'cm^^), we find that the HA plas- 
mon mode has a much higher frequency than the optical 
phonon mode. The experimentally measured plasmon 
dispersion in ref. |l^ is also much higher than the WC 
results shown in Fig. 3 indicating that the WC physics 
is unlikely to be playing a role here. 

In considering the transition to the low density WC 
phase (and in particular whether the plasmon, or more 
appropriately the longitudinal optical phonon, mode in 
the electron crystal phase could be observed in low den- 
sity 2D GaAs based electron systems) we show in Fig. 
4 our best (approximate) estimate for the 2D electron 
liquid/crystal phase diagram (appropriate for 2D GaAs 
based electron systems) in the density-temperature (n, T) 
plot. The classical WC - electron liquid (first order) 
phase transition line is defined by F — F^,, where F — 
{V)/{T), where {V) = 6^(7™) is the classical mean 
potential energy and (T) is the mean classical kinetic en- 
ergy given by 
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FIG. 4. Calculated 2D electron liquid - Wigner crystal 
phase diagram in the density - temperature (n, T) plot ( ap- 
propriate for GaAs systems). The QWC (CWC) represents 
the quantum (classical) Wigner crystal. 

n J {2.TT ) 2m 

where np{p) is the finite temperature Fermi distribution 
function, and Fc is found from numerical (molecular dy- 
namics) simulationsllj to be Fc ~ 120. The high temper- 
ature region above the F = Fc line (dashed line), F < Fc, 
in Fig. 4 is a classical Coulomb liquid (for T > Tp) 
whereas the low T region below the F = Fc line, F > Fc, 
is the classical Wigner crystal (CWC), for T > Tp. At 
zero (low) temperature quantum fluctuations produce a 
(first order) quantum phase transition between the (low 
density) electron crystal phase and the (high density) 
electron liquid phase. At T = this transition occurs 
at = r* where the critical r* is foundtj by quantum 
Monte Carlo simulations to be around r* — 37. The 
high density region (r^ < r*) in Fig. 4 to the right of the 
Ts = r* (at T — 0) line is the quantum liquid whereas the 
low density region (rg > r*) along the T = line is the 
quantum Wigner crystal (QWC). We have produced the 
(n, T) phase diagram in Fig. 4 using a simple mean field 
theory and Fade approximation, which smoothly inter- 
polates between the classical result F = Fc as n ^ and 
the quantum result = r* as T 0. The experimental 
regime explored in ref. |l^ is shown as the shaded block 
in Fig. 4, and we have added the T = Tp line to crudely 
separate the classical and quantum regimes. From Fig. 4 
it is clear that the samples used in ref. are essentially in 
the quantum electron liquid region — increasing T (both 
decreasing n and T) one should be able to explore the 
classical liquid (Wigner crystal) region. We propose that 
plasmon experiments be carried out in the lowest density 
samples at higher temperatures (T > AK) in order to 
test the quantum-classical crossover in the 2D plasmon 
properties. 

Finally in Fig. 5 we show our predicted plasmon 
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FIG. 5. The plasmon spectral weight (or the loss func- 
tion, lm[l/e{q,u})]) at a low 2D electron density, ra — 18.4 
(n = 0.85 X 10^ cm~^) for various temperatures at a given 
wave vector q — lO^cm"^. The Fermi temperature of this 
system is Tp = 350mK. We include an impurity level broad- 
ening of 0.05meV in this calculation. 

behavior in the classical (T > Tp) regime, which should 
be easily achievable (see Fig. 4) in the lowest density 
samples of ref. |ll] simply by raising the temperature to 
3 — 4K. One should observe spectacular thermal level 
broadening of the plasmon peak (as well as a temperature 
induced high energy blue shift of the plasmon energy) as 
shown in Fig. 5. This spectacular level broadening will 
indicate the quantum-classical plasmon crossover — the 
current data in rcf. |l^ are all essentially still in the quan- 
tum regime. Observation of the WC collective modes 
will, however, require samples with substantially lower 
carrier densities (and lower temperature measurements) 
as indicated in our phase diagram in Fig. 4. 

We conclude by summarizing our results: We have ob- 
tained reasonable agreement with the measured plasmon 
dispersionEiJ in low density 2D systems by carrying out 
a complete quantitative calculation of 2D plasma disper- 
sion. We find considerable cancellations among various 
physical mechanisms (e.g., between non-local effects and 
finite width corrections and between thermal effects and 
local field corrections) in the plasmoit-dispersion leading 
to the observed apparent agreemenlllil between experi- 
ment and classical 2D plasma dispersion formula. We 
also study the quantum-classical as well as crystal-liquid 
crossovers in the 2D plasmon behavior and propose spe- 
cific experimental studies to test our theoretical predic- 
tions. 

This work is supported by the U.S.-ONR. 
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